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Abstract 

We calculate fragmentation functions for a b-quark to fragment into 
color-singlet P-wave bound states cb in the Heavy Quark Effective The- 
ory with the exact account of 0(1/ nib) corrections. We demonstrate an 
agreement of the obtained results with the corresponding calculations 
carried out in quantum chromo dynamics. 

Introduction. 

The study of heavy quarkonia properties is a subject of much current interest for under- 
standing of quark-gluon interaction dynamics. _B c -mesons, consisting of b- and c-quarks, 
hold a unique position in the heavy quarkonia physics In the first place, B c -mesons 
consist of two heavy quarks, so the predictions of potential models refer to J/^>- and T 
-mesons as well as to B c -mesons g]. In the second place, B c -mesons are constructed 
from the quarks of different flavors and masses, what essentially determine their decay 
characteristics ||. 

Production of mesons with heavy quarks in e + e~, 77 and pp-interactions may be de- 
scribed in nonrelativistic perturbative quantum chromodynamics. At present, two mecha- 
nisms were investigated for the production of _B c -mesons: the recombination and the frag- 
mentation. In the first case, _B c -mesons are formed from heavy quarks, produced indepen- 
dently in hard subprocess. The fragmentational mechanism demands the pair production 
of b- or c-quarks in hard subprocess with their subsequent fragmentation to _B c -mesons 
(b — ► B c c,c — ► B c b). The relative contributions of these mechanisms in the production 
cross-sections are different in various reactions. In e + e~ -annihilation only quark fragmen- 
tation is essential [Q . In pp 77 -interactions the fragmentational mechanism prevails also 
for _B c -mesons production with large transverse momenta (5[ |6| . In the range of small trans- 
verse momenta the recombination is dominant and determines the total cross section of 
-B c -meson production in 77 and pp interactions. But experimental conditions of i3 c -meson 
discovery are more perspective in the large transverse momenta domain. So, the study of 
heavy quark fragmentation into £? c -mesons attracts considerable interest. An approach for 
the calculation of fragmentation functions Di^ B e , D c ^s c b in nonrelativistic perturbative 
quantum chromodynamics was suggested in |7j . 

At the same time, in the last years there was suggested, based on QCD, the Heavy 
Quark Effective Theory (HQET) j|, |^| for description of heavy hadrons properties. HQET 
make it possible to obtain the finite analytical result with some accuracy even for com- 
plicated processes of quark- gluon interaction. In this approach the matrix elements of 
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different processes may be decomposed on degrees of two small parameters: the strong cou- 
pling constant a s (mq) and Aqcd /rug, where tuq is the mass of the heavy quark. In the 
limit niQ — > oo the effective lagrangian, which describes the strong interactions of heavy 
quarks has an exact spin- flavor symmetry ]8| . HQET is successfully used for investigation 
of exclusive and inclusive hadron decays ||7 Recently it was shown that HQET may 
be used for study of b-quark fragmentation into S-wave pseudoscalar and vector mesons 
and the corresponding nonpolarized fragmentation functions were calculated. The HQET 
calculation of the b-quark fragmentation into the transverse and longitudinal polarized S- 
wave -B*-mesons have been made in [ pT| . In this work we have calculated the fragmentation 
functions of b-quark into P-wave color-singlet states (cb) to the next to leading order in the 
heavy quark mass expansion using the methods of HQET. 



1 Fragmentation functions into P-wave c6-mesons. 

Heavy b-quark may fragment into bound states of two heavy quarks (cb - states) with 
orbital momentum 1=1. There are four such states: 1 Pi, 3 Pj (J=0,l,2). Heavy quark 
fragmentation functions into P-wave B c - mesons were calculated by Chen fl2]| and Yuan fl3|| 
in QCD, but the results of their calculations disagree. Let carry out the similar calculation 
of fragmentation functions in the HQET. Let q — m^v + k is 4-momentum of virtual 
heavy quark, p\ = (1 — r)Mv + p and pi = rMv — p are 4- momenta of b— and c-quarks 
correspondingly; p is 4-momentum of relative motion. Let also I = k — p is 4-momentum 
of the virtual gluon and k is the residual momentum of the fragmenting heavy quark. 
Fragmentation functions for the process b — ► B c + c are determined by the next expression 
§: 



D ^ = TTT^ / ds9 \ 3 lim #FFT2 - (!) 

where M = rrif, + m c is the mass of i? c -meson, T is the matrix element for production 
B c + c from an off-shell 6*-quark with virtuality s = q 2 , and T is the matrix element for 
producing an on-shell b-quark with the same 3-momentum q. The calculation can be greatly 
simplified by using the axial gauge with gauge parameter — (1,0,0, —1) in the frame 
where q^ = (q , 0, 0, ^fq% - s): 



iO 



k a n x + k x n 
9<t\ ; 



(k-nyy ^ 

The part of amplitude T that involves production of the virtual &*-quark can be treated 
as an unknown Dirac spinor T. In the limit q — > oo, the same spinor factor T appears in 
the matrix element To = Tv(q), what leads to cancellation of this factor T in (1). 

Let consider the fragmentation of b-quark into color-singlet bound state (cb) 1 P\. The 
amplitude of such process involves the spinor factor v{pi)u{p2). To project the pair of 
quarks on l P\ bound state we have used the next substitution 



v(P2)u(Pi) -> VM- 75^ ■ (3) 

2m c Zrrib 
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The HQET Lagrangian, including the leading and the l/rUf, terms is given by 



L = h v < iv ■ D 



2rrib 



d(iD) 2 - C 2 (v ■ iD) 2 - 9l gs(JiivG ^ 



h v , 



where 



C x = 1, C 2 = 3 



a s (mb) 



2, C 3 = 



a s (mb) 



' (33-2re ,■) 



(4) 



(5) 



All of these coefficients are equal to 1 at the heavy quark mass scale p = mb- When we 
concerned the fragmentation into P-wave mesons, it is necessary to decompose the projecting 
operator (3) and the gluon propagator (2) on the relative motion momentum p. Using the 
Feynman rules, derived from HQET Lagrangian (4), we may write the full fragmentation 
amplitude into 1 P\ -state Jl4|: 



6m r mh O0n 1 1 



(0) 



\u(p')j u (m c v — p — m c )j 5 (m b v + p + rrib)[v^ + —^-(p + k)^- 
l zrrib 

2mi 2mi V / 



.^. (p + fcK + ^ (p _ fc)A] i^ r 



where e a (L z ) is the polarization vector of Pi -state. To calculate amplitude (6) it is 
convenient to divide it into two parts on the degrees of small parameter 1/m;,. When 
nib — > oo in vertex function and in the propagator of heavy quark, we obtain the main 
contribution to the fragmentation amplitude of b-quark in the form: 

MitfPi) - ^ as2 5!S *Muip')W^ 5 T, (7) 



W a = (s - m 2 b ) 



3r 2 (s — m 2 ) 3 

(V + l)7 a - ^—^rh{v + l)7 a 

n ■ k 





k ■ v 


+ Amk a 


1 + n 


k ■ n 



(v-iy 



—2Mr(s — m?)u Q T n (v ~~ 1) + 2Mr(s — mi)n a - rv7"(^ — !)■ 

?i • k (n ■ k) 

All calculations of the fragmentation functions, which are rather complicated, were done 
by means of the system "REDUCE". Substituting (7) into (1), we obtain: 

D 1 {^P x )[y)=N x { ^^{^-^ + 40^ + 96), N x = "^ffff ■ («) 



where y = (1 — z + rz)/rz is the so called Yaffe-Randall variable Pq[ , and r = m c /M. 
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The amplitude of bound cb state n 3 Pj production may be derived from (6), changing 
75 — > e(S z ), where e fl (S z ) is the spin wave function: 



iM(n 3 Pj) = y/ ^ Ma " R[(0) V < 1,L Z ;1,S Z \J,J Z > e* p (S z )e* a (L z ) (9) 
6m c m\ ) ^ 

^~{^~ 9 ^ v + ^-y)«b')7 ,y ( m ct' - p - m c )-fp(m b v + p + mb)[v^ + ^~(P + k ^~ 

C 2 , . .C 3 „ X/ , s + 

-v 



where we have expressed the Clebsch- Gordon coefficients and e%(S z ), e* (L z ) by the bound 
state polarizations [|f4] : 

-^(9a/3 ~ VaVf}), J = 

J2 <l,L z ;l,S z \J,J z >e* p (S z )e* a (L z ) = { ^e a px p v x e* p {J z ), J=l (fO) 
S*,L* e afi {J z ), J = 2. 

Taking in (9) only the terms of leading order on 1 /mj, and doing necessary differentiation 
on p Q , we have obtained: 



^) = ##E <l,L z ;l,S z \J,J z >e*p(S z y a (L z ) (11) 

S0'){(s - mg)(l - -^n) lal(} - Ak a M(l + -^fi) lp - ArM 2 ( ^—^ J n Q n 7/3 }(l + t))r. 

This amplitude determines the basic contribution to fragmentation functions of b-quark 
into 3 Pj- state: 

Dx(n 3 P )(y) - ^^^(j/ 4 - 4y 3 + Sy 2 + 32), (12) 
ry° 

D 1 (n 3 P i )(y) = ~ 1)2 (3y 4 - Ay 3 + 16y 2 + 48), (13) 

ry° 

D 1 (n 3 P 2 )(y) = Ni 2 ° iv ~ 1)2 (y 4 + Ay 2 + 8). (14) 

The heavy quark fragmentation functions into P-wave mesons were calculated in ]l3[| 
using the QCD. Our results (8), (12)-(14) coincide with the calculations of jf3| in the leading 
order on 1/mj,. Let consider calculation of the next-to-leading order contributions 0(1/ mb) 
for 1 P\ -state. First of all, let observe, that the fragmentation functions, determined by the 
amplitudes (9) and (If), involves factor f/[l + r(y — 1)]. Decompose it on degrees of r, we 
have the correction to (8): 

D 2 (n 1 P 1 )(y) = jVi (1 ~/ )3 (9y 4 - Ay 3 + AOy 2 + 96), (15) 
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Without any calculations we can obtain the propagator correction 0{l/m b ) to (8), which 
arises as a result of following substitution: 



1 + r(-Ci + -C 2 My) 



(16) 



2m h 2m b 

This correction has the next form: 

A^PiXy) - N ^ V ~ y ^ (V - 4 ^ 3 + 4 °y 2 + 96)(-2Ci + C 2 My). (17) 
The vertex correction of necessary order in (6) is calculated in a more tedious way: 

\/AnMa s , d 



iM^X^Pi) = V " ,1V1 " S R' n p(0)e*jL z ) — {u(p')^( mc v - p - m c ) lb (18) 

r r?, 

(mfc{i + p + m b ) 



£t b {p + *>" - £k v ■ {p + k)v » + p) \ 



l+« r i 1 



2 wfcp 



It is natural to perform it as a sum of several terms. Putting p = in the square brackets, 
we obtain: 

C C C 1 

^ty + k) » - *t b v ■ {p + k)v » - ±tM p - - {p - = (i9) 

An addendum (—C2(k-v)v tl /2m b ) in (19) gives the contribution to fragmentation functions 
which differs only by sign from the similar quark propagator correction: 



D^P^iy) = N\— — -p—(9y 4 - 4y 3 + 40y 2 + 96){-C 2 My). 



(20) 



Two other terms in square brackets of (19) give rise the next fragmentation amplitudes: 

, (g ~ ml) 



*^ert.(« P l) = — R np(°) £ a(^) — ~ -5^2 «(P ) 



3m c m b 



2m b r(s — ml) 2 



2k a Mrh- 



k ■ n 



(21) 



( g TTth^)^ 1 (s - 777,^)^ 

-4k a mk + 2n a nMr 2 V n J + rh-/ a K - b kj a (s - ml) 75(1 + fi)r, 



(jfe-n) 



k ■ n 



M&Xn 1 *) = ^C^e^^^ (22) 



3m c m b 



m b r 2 (s — ml) 



■5 



i{p')\-2k a Mrh ^- m ^ - 8k a Mk - 2(s - m 2 b )k a + 2Mr ^ 



k ■ n 



+rh la ( S m ^ + r(s - m 2 b )k la \ l5 {\ + fi)I\ 

A further part of vertex correction 0(l/mb) appears, when we differentiate the expression 
in square brackets of (18) on p a : 



3m c m b 



R' nP (0K(L z ) 



1 



Kp') Ci {^an-^— -7^ +C 3 (k 



2r(s — ml) 2 
1 1 (a -ml) 



(23) 



k ■ n 2M k ■ n 



^nk-f a + 2-f a S j 7 5 (l + 0)r. 



k ■ n 

The contributions of matrix elements (21)-(23) to fragmentation functions have the 
following form: 



D B (n 1 Pi)(y)=4N 1 



(3^-2^ + 6y + 32)Ci, 



D 6 (n 1 P 1 )(y) = w fi ! )y \ -3y 3 - Ay 2 + lAy - 16)C 3 , 



(24) 
(25) 



D^PtXy) = 2Ni — ^ V [(3y 3 + 6y 2 + Ay + 32)d + 2y(3y 2 - 2y - 16)C 3 ] . (26) 

The small component of heavy quark field also leads to the correction of type O(r) in 
function D(y). Substituting the small quark component propagator to (6) instead of heavy 
quark propagator 



i 1 + v 



1 



2m b 



cr^\k > 



1-v 



v ■ k 2 

we obtain the next amplitude of B c - meson production: 

V4ttM2q!. s 



iMpropXn 1 Pi) 



3rn c mb 



-R' nP {0y a (L z ) 



r 2 (s — ml) 



2\3 



(27) 



(28) 



u(p') 



Ak a Mk + 2k a {s-ml) + 



•— —n la + (s - m 2 b )k-/ a 

k ■ n 



7»(i-«)r. 



This amplitude gives the contribution to fragmentation function of the kind: 



D s (n 1 P 1 )(y) = 6N 1 



y° 



(y 5 + 2y 4 + 3y 6 + l6y -16). 



(29) 



The calculation of b-quark fragmentation functions into 3 Pj states was done in a similar 
way. The results of D(n 3 Pj) calculations are represented in the next section. 
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2 Discussion of the results. 



Let analyse b-quark fragmentation functions into P-wave mesons. As was mentioned earlier, 
in the leading order on l/m& the expressions (8), (12)- (14) coincide with QCD calculations 
of Yuan The correction 0(r) to (8) is determined by the sum of terms (15), (17), (20), 
(24)-(26), (29). Setting C\ = C3 = 1, we obtain the fragmentation function D(n 1 Pi) to the 
next-to-leading order on r: 



c^Kp(0)| 2 (2/-1) 2 



547rr 5 M 5 ry s 
r{3y 5 - Sly 4 + 32y 3 H 



Similarly, for the 3 Pj states, we have: 



D( 3 Po)(y) 



a 2 s \K P m 2 (y-^ 

547rr 5 M 5 ry 8 



{9y 4 - 4y 3 + 40y 2 + 96)- 



- 192y + 96) 



(y 4 - Ay 3 + 8y 2 + 32)+ 



(30) 



; (3y 5 - lly 4 + 392y 2 + 192y - 96) 



D( 3 Pi)(y) = 



\Kp(o)\ 2 (y~V 2 



27irr 5 M 5 



ry° 



(3y 4 - Ay 3 + 16y 2 + 48)+ 



D( 6 P 2 )(y) 



r(3y 5 + 5y 4 + 8y 3 + 32y 2 + 96y - 

i0a 2 Kp(0)l 2 {y-i? 



(y 4 + Ay 2 + 8)+ 



15 v y 



27nr 5 M 5 ry 8 
y 4 + 36y 3 - 164y 2 + 240y - 120) 



(31) 



(32) 



(33) 



It follows immediately from (30)- (33), that our results coincide with the calculations of 
Yuan Jl3| with the accuracy 0(r), if we take into account, that the nonperturbative factor 
of (30)-(33) in jl3| contains reduced mass /1 contrary to our factor rM. Using obtained 
fragmentation functions (30)-(33), we may calculate the fragmentation probabilities of cor- 
responding (c6)-mesons 



P b ^- cb {n 1 P 1 )= f dzD b ^{n l P x ){z,^ 
Jo 



') 



(34) 



^Kp(0)I s 
5Airr 5 M 5 



■ In'/ 



(1-r) 



(21 + 32r + 110r 2 - 184r 3 - 387r 4 - 168? 



+ 



210(1 - r) 



y (1032 + 4497r + 21353r 2 - 2762r 3 - 65202r 4 



76199r 5 - 3679r b 



>(n 3 P ) 



a 2 Kp(G)l 



2 . 



54irr 5 M 5 



•In j 



3(1 -r) 



(3 + 8r + 98r 2 + 1496r 3 - 1221r 4 - 960r 5 ) + (35) 
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630(1 -r) 7 
P 5 ^e 6 (n 3 Pi) 



1 



(360 + 1767r - 14977r 2 



180778r 3 + 158658r 4 - 427697r 5 



19849r b 



"?Kp (0)1 5 



277rr 5 M 5 



r In r 



(1 



(3 + 16r + 10r 2 + 112r 3 - 237r 4 - 192r 5 ) + (36) 



210(1 - r) 



Y (348 - 177r + 8419r 2 + 2714r 3 + 14334r 4 - 81769r 5 - 4349r 6 



P b ^ eb (n 3 P 2 



10a 2 K P (0)|^ 
27nr 5 M 5 



r\nr 



15(1 -rf 



(33 - 68r + 130r 2 - 548r 3 - 291r 4 + 24r 5 ) + 



(37) 



-y (1710 + 2697r - 2525r 2 - 10330r 3 - 146760r 4 + 3425r 5 + 583r 6 ) 



3150(1 -r| 

Putting here |i?^(0)| 2 = 0.201 GeV 5 , m c = 1.5 GeV, m b = 4.9 GeV and a s (2m c )=0.38 
(2m c is a minimal energy of exchanged gluon) , we have obtained the numerical value of the 
fragmentation probabilities, which are presented in table. We see, that our integral prob- 
abilities of P-wave cb meson production, founded by means of the b-quark fragmentation 
functions in the HQET are in good agreement with the results of QCD calculations of Jl3||. 



Pb^cbi^Pl) 


Pb^cbi^Po) 


Pb^cb (2 3 -Pl) 


Pb^cbi^P-l) 


6.4 • 10" 5 


2.5 ■ 10~ 5 


7.3 ■ 10" 6 


10.5- 10~ 5 



So the performed calculations show that the Heavy Quark Effective Theory may be 
successfully used for the study of the heavy quark fragmentation. In this approach we may 
systematically take into account the 0(\/mb) corrections in the amplitudes and the proba- 
bilities of the fragmentation, what increases the accuracy of HQET calculations. Moreover, 
it seems more important, that HQET leads to finite analytical answer, when we study com- 
plicated problems in the heavy quark physics. The approach, based on HQET, may be 
used for calculation of heavy quark fragmentation functions into D-wave mesons, and for 
the investigation of B c -meson hadroproduction. 

We are grateful to Faustov R.N., Kiselev V.V., Likhodcd A.K. for useful discussions of 
B c meson physics and the Heavy Quark Effective Theory, and to Braaten E., Cheung K., 
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